Spatial oscillations in the spontaneous emission rate of an atom 

inside a metallic wedge 

H. J. Zhao 1 and M. L. Dufl 

1 School of Physics and Information Science, 
Shanxi Normal University, Linfen O4IOO4, China and 
2 Institute of Theoretical Physics, Chinese Academy of Sciences, 
P. 0. Box 2735, Beijing 100190, China 
(Dated: November 30, 2010) 

Abstract 

A method of images is applied to study the spontaneous emission of an atom inside a metallic 
wedge with an opening angle of tt/N, where N is an arbitrary positive integer. We show the 
method of images gives a rate formula consistent with that from Quantum Electrodynamics. Using 
the method of images, we show the correspondence between the oscillations in the spontaneous 
emission rate and the closed-orbits of emitted photon going away and returning to the atom inside 
the wedge. The closed-orbits can be readily constructed using the method of images and they are 
also extracted from the spontaneous emission rate. 
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I. INTRODUCTION 



The effect of boundaries and environments on spontaneous emission of atoms is an in- 

n 

teresting subject in the physics of radiation with a long history O ur understanding of 
the subject has been greatly improved by the theoretical and experimental developments 
over the years (26]. For example, Hulet et al in an experiment demonstrated clearly that 
the lifetime of atomic state can be greatly modified by the two metal plates [7]. From a 
practical point of view, once the boundaries are given and specified for a system, we would 
like to know how the spontaneous emission rate is modified. In particular, how does the 
spontaneous emission rate depends upon the position and polarization of the emitting atom 



in the system. Quantum Electrodynamics (QED) re 



ates the spontaneous emission rate of 



an atom to the quantized electromagnetic field [8|, |9|. In principle, the calculations of the 
spontaneous emission rate of an atom in a cavity become therefore the search for the quan- 
tized electromagnetic filed of the specified system, the evaluation of the matrix elements 
and the summation of the matrix elements squared. In practice, however, the above process 
can be complicated and difficult to apply even for a simple system such as a metallic wedge 
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11] . Therefore there is a need for alternative methods. 



In this article we will apply a theoretical method of images to study the spontaneous 
emission rate in a metallic wedge. In an interesting application, the method of images 
has recently been used to study the quantum interference at corners[12]. We will show the 
method works for a wedge with an angle ofn/N, where N is an arbitrary positive integer. As 
examples, we have worked out the details for wedges of 7r/3 and 7r/4 corresponding to N = 3 
and N = 4. We find the method of images gives a formula for the spontaneous emission rate 
which is consistent with the one obtained by Rosa et al involving long derivations using QED. 
The method of images is not only easy to apply, it also provide a way to understand the rate 
formula by associating each term with an image and a closed-orbit. For other systems, we 
have shown the large scale oscillations in the spontaneous emission rate [jij, [jjj] are similar to 
;he oscillations for photo-ionization and photo-detachment in the presence of external fields 
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19) . Using scaling variables, we have demonstrated the direct correspondence between 
the oscillations in the spontaneous emission rate and the closed-orbits of photon I3I, 
Here we will use the same method to analyze the oscillations in the spontaneous emission 
rate in a wedge with an angle of 7r/3 and 7r/4 respectively. We will show there are four 
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distinct oscillations for a it/ 3 wedge, and there are six distinct oscillations for a 7r/4 wedge. 

The paper is organized as follows. In the next section, we discuss the construction of 
the images and the calculation of the emission rate. The method of images is general 
in the sense that it can be applied to a wedge with an opening angle of tt/N, where N 
is an arbitrary positive integer. In section III, we take a 7r/3 wedge and a tt/A wedge 
as examples and discuss the oscillation patterns of the spontaneous emission rate as the 
position of the emitting atom inside the wedge is changed. We then analyze and extract the 
oscillations using transformations. We will demonstrate the direct correspondence between 
the oscillations and photon closed-orbits inside the wedges. Conclusions and perspectives 
will be given in section IV. 



II. THE METHOD OF IMAGES AND THE EMISSION RATE 

The cross sectional view of a metallic wedge is shown schematically in Fig.l. The wedge is 
assumed to be made of perfect metal and an atom is placed inside the wedge at the position 
of the black dot. Let 9 be the opening angle of the wedge. Let the z-axis be the wedge 
axis. A cylindrical coordinate system (p, <fr, z) is defined such that the distance between the 
atom and the z-axis axis is p and the declination angle relative to the downside plane is (p. 
A Cartesian coordinate system (x, y, z) is also useful. We define the x-axis in the bisector 
plane of the opening and the y-axis perpendicular to it. For the sake of simplicity, the dipole 
moment of the atom in the following is assumed to be parallel to the wedge axis. 

The method of images has been applied to discuss the spontaneous emission rate of an 

LI 

atom near a mirror and between two parallel plates [9] . The analogy of an radiating antenna 
to an atomic emission process is explored in the method of images. For an emitting atom 
near a plane mirror, it has been demonstrated that the spontaneous emission rate can be 
inferred from the power radiation rate of an antenna |9] . Consider the radiation damping of 
a dipole antenna with the dipole moment cfe~ lWo< , where ojq is the frequency. The radiation 
damping rate is defined as Wd = P/U, where P is the average radiation power of the antenna 
and U is the mechanical energy of the antenna. Let f denote a vector of a point relative to 
the dipole position, the radiation damping rate can be written as [9| 

W d = ^jMd* ■ E(f= 0)] = ^{Im[d* ■ E (f = 0)] + Im[ef ■ E ref (f= 0)]}, (1) 
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where the electric field E at the position of dipole antenna has been decomposed into a 
direct part E and a reflected part E re f. The direct electric field is 

* = x <*> x %> + |d - 3f " • d)l x W - wy^"^' (2) 

where c? is a unit vector in the dipole direction and Eq is the vacuum permittivity. If the 
dipole antenna is at distance 1/2 from a plane metal boundary and is parallel to the metal 
surface, the effect of the conducting boundary can be replaced by an image dipole and the 
reflected electric field E re f of the dipole can be calculated directly from the imaging dipole. 
The imaging dipole direction is opposite to the original dipole direction and is located at a 
distance 1/2 behind the metal surface. The reflected field at the position of the antenna is 

where / is the distance between the image and the original antenna. 

When a dipole antenna is placed in a metallic wedge as shown in Fig.l, the reflected field 
can be calculated by using a set of images and the reflected filed is a sum of contributions 
from each image of the set. We find when the opening angle of the wedge is 9 — tt/N, where 
TV is an arbitrary positive integer, the number of the images is finite. The images can be 
constructed using a technique derived from a wave-packet evolution study inside a wedge by 
Robinett jsO | - If the emitting atom's cylindrical coordination is at (p,<p,z), the first image 
with a dipole —d will appear at (p, —<f>,z). The original source and the first image can be 
rotated by 29 clockwise to give the second image located at (p, (f> — 29, z) with a dipole d and 
the third image at (p, —0 — 29, z) with a dipole —d. The original source and the first image 
are rotated by 49 clockwise to give image 4 located at (p, <p — 49, z) with a dipole d and image 
5 at (p, —(f) — 49, z) with a dipole —d. This process can be continued until the last two images 
are obtained by rotating the original source and the first image by 2(N — 1)9 clockwise. For 
9 = tt/N, where iV is an arbitrary positive integer, this rotation process closes on itself. In 
the end, we have one original source and (2iV — 1) images. (N — 1) of those images with the 
same dipole d as the original's are located at (p, <fi — 2nn /N, z), where n = 1,2,..., (N — 1). 
The other N images are located at (p, — <j) + 2nir/N, z), where n — 0,1, ... , (N — 1), and they 
have the opposite dipole — d as the original's. Fig. 2 shows the five images and the original 
source for a 9 = n/3 wedge. 

The reflected electric field at the position of the source antenna in the wedge with an 
opening angle 9 = n/N can be obtained from the electric fields of all the images using Eq.(3). 



The distance between the source at (p,<f>,z) and the image at (p,—</>,z) is /_ = 2psin(0); 
the distance between the source at (p,<fi,z) and the images at (p, (f> — 2nn/N,z), where 
n = 1, 2, . . . , (N — 1), are l n+ = 2psm(nn/N); the distance between the source at (p, 0, z) 
and (p, —(f> + 2nir/N, z), where n — 1, 2, . . . , (N — 1), are /„_ = 2psin(n7r/iV — 0). 
The reflected field at the position of the original dipole antenna can be written as 

E t = -—d\— + — —} e W— <*>t) 

r6/ 4:lT€o kl- (kl-) 2 (kl-)* 1 

CLK c%r^. r 1 2 1 ^ o i(kl n --w t) 



dj^tij^ — " |( 



47re ^ kin- (kin-) 2 (kl n -) 3 
1 7 1 

_ j _|_ __: ] e »(*w-^o*)j _ (4) 

/c/ n+ (kl n+ ) 2 (kl n+ ) 3 

Inserting Eq. (TJJ and Eq. (J2]) into Eq. ([[]) and after some manipulations, we get the damping 
rate inside the wedge 

o o N-i 

W d = W {1 - -G||[2A;psin0] - - £[G||[2fcpsin(^ - 0)] - G\\[2kp sin(^)]]}, (5) 

n=l 

where we have defined 



_ sinS* cosS sinS* 

G l{ (S) = — + — (6) 



and Wo = | 4^ w is the damping rate in free space [9|. 

The spontaneous emission rate in a wedge 7 can be inferred from Eq.(5) as 

7 = 7o{l - fapkpsmct)} - | ^[G^kpsmQ - 0)] - G^kpsinQ)]]}, (7) 

n=l 

where 70 is the spontaneous emission rate without the wedge. It should be noted that the 
method of images has its limitation and the correct value of 70 from QED will have to be 
used in Eq.(7)|g, bl|- 

Rosa et al studied the emitted power of an atom in a perfect wedge with an arbitrary 
opening angle using QED. The emitted power is proportional to the spontaneous emission 
rate. The formula derived by Rosa et al using QED involves special functions and appears 
quite complicated. But when the opening angle of the wedge is restricted to 9 = n/N, where 
N is an arbitrary positive integer, the QED formula can be simplified considerably (see Eq. 
(46) in Ref. [10]). We find Eq.(7) using the method of images is consistent with the result 
of Rosa et al using QED for a 9 = tt/N wedge. It is not easy to see why the complicated 
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QED formula involving special functions obtained by Rosa et al can be reduced to a simple 
form similar to Eq.(7) when the opening angle of the wedge is 9 = tt/N, where N is an 
arbitrary positive integer. The method of images seems to suggest that the origin of the 
simplification is in the symmetries of the finite number of images. When the opening angle 
is tt/N, where N is an arbitrary positive integer, the 2N — 1 images and the original source 
form a set, this set is invariant when it is rotated by a multiple of 2n/N either clockwise 
or anti-clockwise. Furthermore, the set can be changed to a another set by an inversion 
through either the upper plane or the lower plane of the wedge. The effect of the inversion 
is equivalent to reversing the dipole direction of the original set. However, the rate formula 
in Eq.(7) is invariant under the above operations as is apparent from the derivation of the 
method of images. We believe the reduction of the QED formulas to the simple form should 
be related to the above mentioned symmetries. 



III. OSCILLATIONS IN THE EMISSION RATE AND CLOSED-ORBITS 



We now analyze the spontaneous emission rate of an atom inside a tt/N wedge. In Fig. 3 
we show the pattern of the spontaneous emission rate described by Eq.(7) for a 9 = tt/3 
wedge. The overall spontaneous emission rate pattern resembles "an egg tray" turned upside 
down. There are two main features in the pattern in Fig.3. They are the suppression of the 
emission rate near the wedge surface and the regular oscillations of the rate inside the wedge. 
For the suppression of the rate, we note in the small S limit, G\\(S) can be approximated by 
(| — j^S 2 ). Therefore the rate goes to zero quadratically as the emitting atom approaches 
the wedge surface. The quadratical nature of the rate suppression can be seen in Fig.4 
in which the emission rate is shown along two straight lines parallel to x-axis and y-axis 
respectively. 

As we will explain next, the oscillations are related to some special photon closed-orbits 
along which the emitted photon travels away from and returns to the emitting atom inside 
the wedge. The oscillations in the cross sections for photo-ionization and photo-detachment 



in the presence of external fields are success 



going away and returning to the nucleus [15 



ully related to the closed-orbits of photo-electron 



- Il9j . Recently the oscillations in the spontaneous 



emission rate have also been related to photon closed-orbits going away from and returning 
the emitting atom in two systems l^, Jjj]. In the following, we will first analyze the oscilla- 
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tions in the emission rate for a 7r/3 wedge and will demonstrate the correspondence between 
the oscillations and some photon closed-orbits. A similar analysis is then performed for a 
7r/4 wedge which confirms our understanding. 

we define a scaling variable a by p = ctAo, where Ao = 27r//c is the wavelength in vacuum. 
We now examine the spontaneous emission rate in the wedge as a function of a with an 
arbitrary but fixed of the emitting atom. Define the modified Fourier transform (FT) 

7(5) = / [ 7 (a) - l0 }ae iaS da, (8) 

where 7(a) represents the dependence of the emission rate as a function of the scaling 
variable a. In our numerical calculations, we took ai = 1.0, «2 = 26.0 and the integration 
step size Act = 0.01. As an example, for <ft = 7r/9, the emission rate as a function of the 
scaling variable a for the 9 = 7r/3 wedge is shown in Fig. 5 (a). One can see the rate in the 
wedge oscillates around the free space rate and the oscillation amplitude is reduced as a 
increases or the emitting atom is further away from the wedge axis. The absolute value of 
the calculated modified Fourier transform is presented in Fig.5(b). There are clearly four 
large peaks in Fig. 5(b). These four peaks imply the emission rate in Fig. 5 (a) contains four 
main frequencies. We now provide an explanation of the physical meaning for the peaks and 
give a quantitative description for the peak positions and peak heights. 

We note when S is larger than about it, G\\(S) = is an accurate approximation of 
G\\(S). When this approximation is made in Eq. (7), the emission rate inside the wedge 
takes the following form 

7(a) = 7o + 5^(^^)sin(S'°a + ^), (9) 



j 
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which is in the form of closed-orbit theory for absorption spectra |15l-ll8l|. Eq.(9) can be 
interpreted using closed-orbit theory. 70 is a background rate associated with the initial 
emission process in which the emitted photon leaves the atom and never returns. It is 
equal to the spontaneous emission rate in free space. The sum in Eq.(9) is over all closed- 
orbits of emitted photon going away from and returning back to the emitting atom; the 
emitted photon travels in straight paths in the wedge, and it follows the laws of reflection 
when hitting the metal boundary of the wedge. S® = hi® is the "action" around the jth 
closed-orbit, k is the wave number, and 1® is the geometric length for the jth closed-orbit 
corresponding to the case a = 1 or the distance between the emitting atom and the wedge 



axis is just equal to one wavelength; since each reflection by the wedge surface produces a n 
phase shift, the phase ipj = —rrijiT where m,j counts the number of reflections by the wedge 
surface; the amplitude (j^t^) provides a measure of the intensity of the returning group of 
photon wave, and in the present case, it is inversely proportional to the length of the orbit. 
The peak positions in the Fourier transform of Eq.(8) are therefore given by the action S® 
of photon closed-orbits with corresponding peak heights proportional to gj/l®, where gj is 
a degeneracy factor counting closed-orbits with identical contribution to the peak. Both S® 
and 1® are evaluated for the system size corresponding to a = 1 or p = X . 

For a general wedge, we must launch trajectories in all direction from the location of 
the emitting atom at (p, 0, z) to search for all closed-orbits which leave from and return to 
the emitting atom. Each trajectory follows a straight line until it reaches the wedge and is 
reflected. When the trajectory returns back to the location of emitting atom, a closed-orbit 
is then formed. The closed-orbits must be on the cross-section plane of the atom because 
of the propagation rules of trajectories. If a trajectory leaves the atom off the cross-section 
plane, it can not return to the atom to form a closed-orbit. For a wedge with an arbitrary 
opening angle, the number of closed-orbit can be large. But for a n/N wedge, we can find 
all the closed-orbits from the images constructed earlier. We find there is a one to one 
correspondence between the images and the closed-orbits. 

In Fig.6, we show this correspondence for a 9 = n/3 wedge. There are five images labeled 
from 1 to 5 clockwise. There are also five closed-orbits. The first segment of the closed- 
orbit is given by the straight line connecting the corresponding image and the source. The 
geometrical length of the closed-orbit is equal to the distance between the atom and the 
image. The first closed-orbit leaves the emitting atom downward, it returns to the atom 
after being reflected back by the lower wedge surface; the second closed-orbit is first reflected 
by the lower wedge surface and then by the upper wedge surface before it returns to the 
emitting atom; the third closed-orbit is first reflected by the lower wedge surface, it then 
travels perpendicular to the upper wedge surface, finally it retraces back to the emitting 
atom after being reflected by the upper wedge surface. Similar descriptions can be given 
to the fourth closed-orbit and the fifth closed-orbit. We also find closed-orbit two and four 
have the same shape but their propagation directions are opposite. 

The scaled "actions" of the five closed-orbits depend only on the emitting atom's cylin- 
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TABLE I: The four peak positions and heights numerically extracted using FT are compared with 
the analytic predications discussed in the text for a ir/3 wedge. The analytic results are represented 
by the crosses in Fig. 5(b). Peaks are labeled from left to right with increasing value of the scaled 

action. 

position height 



peak label 


numerical 


analytic 


numerical 


analytic 




4.30 


4.30 


4.37 


4.37 


2nd 


8.08 


8.08 


2.30 


2.32 


3 rd 


10.89 


10.88 


3.42 


3.45 


4th 


12.40 


12.38 


1.50 


1.52 



drical angle and they are 

= 2£;A o sin(0) = 4vrsin(0), 

Si = 2/cA sin(^) = 47rsin(^), 

S° = 2A;A sin(^ + 0) = 47rsin(^ + 0), 

o0 _ c0 

S% = 2k\ sin(^ - 0) = 4tt sin(^ - 0). (10) 

Since each reflection by the wedge induces a phase loss of ir, the phases corresponding to 
the five closed-orbits above are 

fl = ¥5 = -7T, 
<£> 2 = (f 4 = -27T, 

(p 3 = -37T. (11) 

We have used the formula hgi/lf to compute the i th peak height located at S® for the 
7r/3 wedge. The results from the analytic formulas with = 7r/9 are marked as crosses in 
Fig. 5(b). The degeneracy factor gi is equal to 1 for the first, second and fourth peak because 
one closed-orbit corresponds to each of the three peaks; for the third peak near 10.89, 
because there are two closed-orbits making the same contributions in both the emission rate 
and the FT, the corresponding degeneracy factor is 2. For the purpose of comparing with the 
numerical results, the constant h was determined by matching with the first peak extracted 



TABLE II: The six peak positions and heights numerically extracted using FT are compared with 
the analytic predications discussed in the text for a it/ 4 wedge setting <p = n/12. The analytic 
results are also represented by the crosses in Fig. 8(b). Peaks are labeled from left to right with 
increasing value of the scaled action. 

position height 



peak label 


numerical 


analytic 


numerical 


analytic 




3.24 


3.25 


5.82 


5.79 


2nd 


6.30 


6.28 


3.05 


3.00 


3 rd 


8.88 


8.89 


4.37 


4.24 


4th 


10.85 


10.88 


1.67 


1.73 


5 th 


12.13 


12.14 


1.53 


1.55 


6 th 


12.56 


12.57 


1.45 


1.50 



from FT. We can see that all the peak positions and heights are accurately described by the 
analytic formulas. In Fig.5(b), we also show the closed-orbits correlated with the peaks. A 
detailed comparison is made in Table I between the numerically extracted peak positions 
and heights and the analytic formulas. We find they are in good agreement. 

To confirm our understanding, we also studied the emission rate in a tt/4 wedge setting 
= 7r/12. The results are summarized and presented in Fig. 7 and Fig.8. The images are 
constructed in Fig. 7(a). The pattern of the emission rate inside the wedge calculated using 
Eq.(7) with A^ = 4 is shown in Fig. 7(b). This pattern also resembles "an egg tray" turned 
upside down. The seven images and the seven corresponding closed-orbits are illustrated in 
Fig. 7(c). The emission rate as a function of the scaling variable a is shown in Fig. 8(a). The 
numerically calculated FT is displayed in Fig. 8(b). We find there are six main peaks in the 
FT. The crosses near the peaks mark the results of the analytic formulas for the six peak 
positions and peak heights. The closed-orbits corresponding to the peaks are illustrated by 
the insert figures. From left to right, each and every peak is associated with one closed- 
orbit except the third one. The third peak corresponds to two closed-orbits. For the it/ '4 
wedge, the numerically extracted peak positions and peak heights are listed together with 
the results of the analytic formulas in Table II. Again we find the numerical results and the 
analytic formulas agree well. 



10 



IV. CONCLUSIONS 



In conclusion, we have studied the spontaneous emission of an atom in a n/N wedge 
using a method of images. The method of images is particular simple and the results for the 
emission rate are consistent with the results of QED. From the method of images, we can 
see that the rotation and inversion symmetries should play important roles in the reduction 
of complicated QED rate formulas to simple forms when the wedge opening angle is ti/N, 
where N is an arbitrary positive integer. For a wedg e with an opening angle, which is not 
equal to tt/N, the closed-orbit theory 0, Q, Q, I Til ] suggests the emission rate can still be 
expressed in a form similar to equation (7) although the number of oscillatory terms can be 
large, and the corresponding closed-orbits will have to be searched numerically \v\ . 

We analyzed in detail the oscillations in the emission rate using scaling variables and 
Fourier transformation for the N = 3 and N = A examples. We found only four 

oscillations in a tt/3 wedge and only six oscillations in a 7r/4 wedge. We found a direct cor- 
respondence between the oscillations in the emission rate and closed-orbits. We interpreted 
the oscillations as interferences between the outgoing and incoming photon waves traveling 
along closed-orbits from the atom to the atom. The closed-orbits can be constructed readily 
from the method of images. The method of images provides a complementary perspective 
on the spontaneous emission process for an atom in a wedge, keeping in mind the difference 



between the method of images and the QED method 21]. It is interesting to note that in an 
recent experiment, Fleet et al has put a classical dipole radiating at wavelengths of the or- 



der of a few centimeters inside a wedge-shaped cavity and measured the radiated power 22] . 
They observed that the classical dipole also exhibit effects similar to the Purcell effects. 
However, the classical dipole antenna is not a point source and can not yet simulate the 
emission of an atom in a wedge. On the other hand, the recent observation of the signature 
of photon periodic orbits in the spontaneous emission spectra in laterally confined vertically 



emitted cavities 23j] provides a very encouraging sign that an observation of the photon 



closed-orbits in a wedge is possible in the near future. 
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FIG. 1: Schematic diagram for an emitting atom (black dot) in a wedge with an opening angle 
6 = it/3. Both cylindrical and Cartesian coordinate systems are used and illustrated. The z-axis 
is in the direction perpendicular to the paper. The x-axis is in the bisector plane of the wedge. 
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FIG. 2: The emitting dipole and their images for the = ir/3 wedge. The images represented by 
have dipoles identical to the dipole of emitting atom while the images represented by <g> have 
opposite dipoles. 
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FIG. 3: The pattern of the spontaneous emission rate is formed as the position of the emitting 
atom is varied inside a wedge with an opening angle = ir/3. 70 is the emission rate without the 
wedge. The pattern is like "an egg tray" upside down. 
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FIG. 4: The spontaneous emission rate for the emitting atom at (a) kx = 15 and along the y-axis 
(b) y = and along the x-axis. Dotted lines correspond to Eq.(7) with G\\ from Eq.(6) while solid 
lines correspond to Eq.(7) with G\\ approximated by 
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FIG. 5: (a) The spontaneous emission rate for the = tt/3 wedge as a function of the scaling 
parameter a = p/\o with <p = ir/9 for the emitting atom; (b) FT (absolute value) of the above 
emission rate as defined in Eq.(8). The crosses near the peaks mark the peak positions and heights 
obtained from the analytic formulas discussed in the text. The closed-orbits corresponding to the 
peaks are also displayed and illustrated with the same colors (online). 
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FIG. 6: The closed-orbits and their corresponding images in a 7r/3 wedge. Each closed-orbit 
corresponds to one image. The first segment of the closed-orbit is given by the straight line 
connecting the corresponding image and the source. 
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FIG. 7: The images, the emission rate and the closed-orbits for a 7r/4 wedge, (a) The images 
are constructed according to the method discussed in the text; (b) the pattern of the emission 
rate inside the wedge as a function of the position of the emitting atom; (c) the closed-orbits are 
constructed from their corresponding images. The first segment of each closed-orbit is given by 
the straight line connecting the corresponding image and the source. 
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FIG. 8: Similar to Fig. 5 but for the 7r/4 wedge with (j> = n/12. 
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